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MAGNETIC FIELDS — 

EARTH AND EXTRATERRESTRIAL 

1. INTRODUCTION 

Natural magnetic fields as well as induced fields are important considerations in space 
vehicle design because they can adversely affect onboard equipment and scientific 
instrumentation and can change an orbiting spacecraft’s drag, attitude, and direction of 
motion. For example, the magnetic torques resulting from the interaction between the 
ambient magnetic field and the spacecraft’s magnetic properties may degrade the 
performance of the attitude control system (ref. 1). 

This monograph presents the information on the strength and direction of the natural 
magnetic fields needed for vehicle design. These fields include the geomagnetic field 
(confined to the volume of space about Earth called the magnetosphere), the interplanetary 
field (Sun’s field), and the fields of the planets and their natural satellites. Knowledge of 
these natural magnetic fields is increasing rapidly, but the present state of knowledge limits 
the availability of a definitive analytical model of the Earth’s main (geomagnetic) field to 
about 6.6 Earth radii (synchronous altitude) as presented herein. For the regions beyond 
this altitude, the available information has been reviewed to determine both quantitative 
and qualitative data useful for design purposes. 

Induced magnetic fields encountered during manufacture, testing, and transportation of 
spacecraft can cause undesirable magnetic properties to be introduced into spacecraft. These 
fields warrant careful consideration and control in order to achieve “magnetic cleanliness” 
of spacecraft- the subject of a separate monograph. References 2 and 3 provide some 
information on this subject. 


2. STATE OF THE ART 

Extensive observations of the magnetic field on and above the Earth’s surface have been 
made for more than a century. Recent data from space exploration have changed the 
concept of the geomagnetic field from one of a simple dipolar field to that of 
the magnetosphere illustrated in figure 1 (ref. 4). This figure shows that the solar plasma 
flow (solar wind) compresses the Earth’s main field on the subsolar side and extends it as a 
long tail (magnetotail) in the antisolar region. The interaction of the solar wind with the 
Earth’s main field creates a shock front on the subsolar side. Between this shock front and a 
boundary layer (magnetopause) is a region of solar wind and interplanetary field turbulence 
called the magnetosheath. The interplanetary field is deflected around Earth’s main field at 
the magnetopause which delineates the extent of the Earth’s main field into space. 

Figure 2 (ref. 5) illustrates a useful concept of the solar wind interaction with Earth’s main 
field advanced by V.C.A. Ferraro and S. Chapman in 1931. The solar wind, which is an 
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Figure I . The magnetosphere (ref. 4) 
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Figure 2.-Ferraro-Chapman model of rnagnetospheric compression (ref. 5). 


electrically neutral plasma consisting mostly of protons and electrons, impinges on the 
magnetopause. The interaction of these charged particles with the Earth’s main field causes 
the particles to be deflected away from the magnetosphere, an action which generates a 
current at the magnetopause. The deflected particles are swept around the magnetosphere 
by the solar wind. 

2.1 Earth’s Main Field 

The Earth’s main field is an internally generated magnetic field that resembles the magnetic 
field of a dipole and extends throughout the magnetosphere. The axis of the dipole passes 
through northern Greenland at 78.5° N and 69.0° W and has an inclination of about 1 1.5° 
with the Earth’s axis of rotation (ref. 6). 


2.1.1 Surface Geomagnetic Field 

The surface geomagnetic field consists of the main field generated within the Earth and all 
external field components that are measurable at the Earth’s surface. The surface intensity 
ranges from about 0.3 gauss (3 X 10' 5 T)* at the equator to about 0.7 gauss (7 X 10 s T) at 
the poles as illustrated in figure 3 (ref. 7). Departures of the field from a dipole field are 
called anomalies. These anomalies are both regional and local. Regional anomalies, such as 
the South Atlantic anomally, cause irregularities in the surface field over an area of several 
thousand square kilometers. Local anomalies cause irregularities of about .01 gauss ( 10 3 nT) 
or less over an area of a few square kilometers to several hundred square kilometers. 

Measurements of the geomagnetic field intensity and direction have been made by 
observatories, magnetic survey ships, and more recently by instrumented aircraft in data 
collection efforts. The most recent charts of geomagnetic field data at the Earth’s surface 
for various geographical regions are listed in “Catalog of Nautical Charts and Publications, 
1968” (ref. 8) and are available from the U. S. Naval Oceanographic Office**. Magnetic 
measurements are referenced to geographical or geomagnetic dipole coordinate systems. 
Figure 4 (ref. 9) illustrates the relationship of coordinate systems for the Earth’s surface 
which are based on the average dipole axis (geomagnetic dipole coordinates) and the spin 
axis (geographic or geodetic coordinates).*** 

Magnetic measurements indicate the existence of secular and temporal variations in the 
surface geomagnetic field. 


* Dimensions, units, and conversion factors are defined in appendix B. 

** U.S. Naval Oceanographic Office, Suitland, Maryland 20390. Hydrographic Surveys Department, 
Magnetics Division. 

*** Appendix A presents a summary of the coordinate systems commonly used in magnetic field studies 
and mentioned in this monograph. 
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2.1.2 Secular Variations 


Secular variations are long-term, small magnitude field changes that are measured in terms 
of years. The predictability of these variations is limited because an acceptable theoretical 
model of the sources of the main field has not been developed. The change (growth or 
decay) in the intensity of a secular variation is as large as 150 7 (150 nT) per year (ref. 6). 
Although this change is nominal in a year’s time, the accumulation of several yearly secular 
changes completely outdates magnetic measurements for a particular epoch, i.e. reference 
date. Secular variation charts that indicate field changes as a function of time are issued by 
the U.S. Naval Oceanographic Office. 


2.1.3 Temporal Variations 

Temporal variations are short duration disturbances in the geomagnetic field resulting from 
solar activity and changing relative positions of the Sun and Earth. These variations typically 
have durations ranging from a few seconds to several days and amplitudes from a few 
hundredths to several hundred gammas. 

The frequency of occurrence is determined by the rotation rates of the Earth and Sun. The 
Sun’s rotation can place an active solar area opposite the Earth once every 27 days. These 
active solar areas may persist for as long as 8 to 10 solar rotations. A semiannual increase in 
the magnitude of temporal variations occurs during the March-April and September-October 
periods when the Earth is near the equinoxes. 

The intensity of temporal variations is related to the number of sunspots. The number of 
sunspots in turn changes with the 1 1-year solar cycle as shown in the prediction in figure 5. 
Intensity also varies with the type of temporal variation; diurnal, magnetic storm, or sudden 
impulse disturbance. 

Diurnal variations are smoothly changing day-to-night variations in geomagnetic field 
intensity. The change in any field component caused by diurnal variations is generally less 
than 200 7 (200 nT) (ref. 10). The characteristics of these disturbances are discussed in 
section 2. 1 .4. 

Magnetic storms are intervals of magnetic activity that occur almost simultaneously on a 
worldwide basis. Magnetic storms progress through the characteristic phases illustrated in 
figure 6 (ref. 3). Occasionally the initial phase is preceded by a preliminary reverse phase 
which has a duration of a few minutes. The initial or sudden commencement phase has a 
duration of about one hour during which the total field intensity increases by about 50y (50 
nT). The main phase has a duration of a few hours during which the total field intensity 
decreases by more than 40Oy (400 nT). During the main phase, the horizontal field 
component has a marked decrease in intensity. The greatest change occurs in the polar 
regions and can exceed 1000y (1000 nT). The initial recovery phase has a duration from 
about six hours to about two days during which time the total field intensity rapidly returns 
to near pre-storm levels. 
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Figure 5.— Predicted sunspot numbers. 
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Figure 6. Magnetic storm effect on main geomagnetic field intensity (observation 
of a large magnetic storm in Hawaii, February, 1958) (ref. 3). 


The change in the surface geomagnetic field resulting from magnetic storms can be as large 
as 400 7 (400 nT). The final recovery phase has a duration of several days during which time 
the total field intensity returns to pre-storm levels. The maximum change in the surface 
geomagnetic field is about 100 to 150 7 (100 to 150 nT). A typical surface field change is 
about 40 7 (40 nT) (ref. 1 0). 

Sudden impulse disturbances are short duration, low amplitude fluctuations of the 
geomagnetic field on a worldwide basis as shown in figure 7 (ref. 5). A sudden impulse 
disturbance has a duration of about one hour or less. Amplitudes as large as 
±187 (±18 nT) have been observed for these disturbances at the Earth’s surface. 
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Figure 7. Sudden impulse disturbance (ref. 5). 

Very small magnetic field intensity changes which are oscillatory and of short duration are 
superimposed on the temporal variations described above. These changes are called 
micropulsations and have periods of less than two or three minutes and amplitudes of less 
than a few gamma. Micropulsations with a period less than seven seconds and an amplitude 
less than one 7 (1 nT) have been observed (ref. 1 1). 


2.1.4 Ionospheric Currents and Related Disturbances 

Ionospheric currents result from systems of charged particles moving between 50 kilometers 
and the outer reaches of the ionosphere (approximately 600 kilometers) (ref. 9). Changes in 
ionospheric currents are the main sources of the diurnal variations of the geomagnetic field. 
Because the intensity of the magnetic field resulting from these currents decreases rapidly 
with increasing distance, the ionospheric currents do not make significant contributions to 
the geomagnetic field at distances greater than a few Farth radii. R e , from the Earth’s 
surface (ref. 12). Those currents which flow perpendicularly to the geomagnetic field lines 
are most intense near 1 00 kilometers and those flowing parallel to field lines reach a 
maximum intensity at approximately 150 kilometers (ref. 1 1 I. 


2. 1.4.1 Auroral Currents 

Auroral currents are the polar electrojets and other transient currents associated with aurora 
and changes in the geomagnetic field. They occur mainly in the polar ionospheric regions 
above 67° geomagnetic latitude or about 23° from the geomagnetic poles (ref. 9). However, 
aurora have occasionally occurred in the middle latitudes. Auroral activity has a periodicity 
that is related to the surface activity of the Sun and its 27-day rotation period. According to 
reference 1 1, maximum activity occurs during the spring and fall when the Earth is at the 
equinoxes. The polar electrojet usually flows westward and encircles the polar region with 
part of the electrojet crossing the polar cap and rejoining the circular electrojet in the 
morning sector. Intensity changes in the polar electrojets can cause intense negative 
geomagnetic bays (polar magnetic substorms whose recorded curve resembles the shape of a 
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geographical bay in a coastline). Intense negative excursions as great as 2000 7 (2000 nT) 
have been observed. These excursions as shown in figure 8 (ref. 11) are typically about 
1500 7 (1500 nT) and have a duration of about 30 to 120 minutes. An intense negative bay 
starts in the midnight sector of the electrojet and propagates in all directions. 



A. Observation at Sitka, Alaska 



B. Observation at College, Alaska 


Figure 8.- Horizontal component of surface geomagnetic field during a polar magnetic 
substorm (ref. 11). 



2. 1.4. 2 Equatorial Electrojets 


The equatorial electrojets are laterally limited, relatively intense currents located in the 
ionospheric region above the magnetic equatorial zone at an altitude of about 95 to 110 
kilometers. They generally flow in a west-to-east direction in the sunlit hemisphere. 
Observations made during rocket flights indicate the presence of lower and upper current 
layers. Heppner in chapter V-2 of reference 1 1 cites Cahill's interpretation of the rocket 
data as indicating the presence of a west-to-east electrojet about 20 to 25 kilometers above 
the lower current layer. In the same chapter Heppner cites subsequent research and analysis 
by Maynard which led him to substanially the same conclusion. 


2 . 1 .4.3 Solar Currents 

Solar currents in the ionospheric E region are called solar quiet or quiet day solar daily 
variation (S q ) currents. The daily magnetic variation induced by changes in the S q currents 
is described in terms of an average of the 5 days per month that had the least magnetic 
variation (quiet days) or more commonly in terms of an average of quiet days selected from 
each month. S q current systems are formed in the subsolar and antisolar hemispheres. These 
systems are stationary with respect to the Sun but have a 24-hour period with respect to an 
observer on the Earth. Increased activity on the Sun’s surface causes an increase in the ion 
density and conductivity of the ionosphere, thus increasing the intensity of the quiet-day 
variation currents in the subsolar hemisphere. The S q currents in the middle latitude regions 
cause variations in the surface geomagnetic field of about 20 to 40 7 (20 to 40 nT). The S q 
currents located near the magnetic equator (Huancayo, Peru) cause surface variations of 
about 100 to 200 7 (100 to 200 nT) (ref. 10). 


2.1 .4.4 Lunar Currents 

Lunar (L) currents in the ionospheric E region are similar to quiet-day variation currents but 
are only about 1/30 as intense. L currents are generated by the gravitational attraction of 
the Moon, causing movement of ion masses in the ionosphere. These movements cause 
variations in the surface geomagnetic field of about 17(1 nT) with a period of 12 hours. 


2.1.5 Ring Currents 

Ring currents are produced by the movement of electrons and protons in circular spatial 
distributions within the magnetosphere. The intensity and spatial distribution of ring 
currents are not well defined and therefore these sources of geomagnetic field variations are 
not usually included in field calculations (ref. 12 ). 


2.1.6 Magnetic Field at Synchronous Altitude 

The equatorial magnetic field at the geocentric distance of about 6.65 R e consists of the 
steady-state geomagnetic field plus the temporal variations of this field. On magnetically 
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quiet days, diurnal magnetic field disturbances due to external field sources cause a 10 to 
15% variation in the total magnetic intensity. During magnetic storms, the field at 
synchronous altitude can increase by 50%, decrease to zero, and even reverse in direction. 

Figure 9 shows the intensity of the field components as a function of longitude at 
synchronous altitude. (These curves are derived from a 24 -term spherical harmonic 
expansion and reflect only that portion of the field due to internal sources.)* The curves are 
reasonably accurate for magnetically quiet-day conditions along either the dawn or dusk 
meridian; the values at these meridians differ by about 1 5%. 

Figure 10 (ref. 13) illustrates the diurnal variation of the field components at the longitude 
of Hawaii on a magnetically quiet day. The vertical and east component showed little 
variation during the day. The horizontal component has a variation of about 
± 15 7 (± 15 nT). The magnetic field intensity is maximum at local noon (2200 hours 
Greenwich mean time at Hawaii) because of the steady -state compression of the 
magnetosphere by the solar wind and at minimum intensity at local midnight because of 
expansion of magnetosphere in the antisolar direction. 

Figure 11 (ref. 13) presents field measurements at the longitude of Hawaii made by the 
ATS-1 synchronous satellite during a magnetic storm on 13, 14 and 15 January 1967. 
Shortly after local noon on 14 January the magnetopause dipped below synchronous 
altitude. This event is indicated by the discontinuity in ATS-1 measurements-marked 
“boundary crossings” on the figure. During a two-hour period the horizontal component of 
the field reversed in direction with a maximum magnitude of 140 7 ( 140 nT) and. 
within an hour, reached a peak of +180 7 (+180 nT). Shortly after local noon of 
15 January, magnetic conditions had returned to the typical quiet-day variations. 

2.1.7 Cusp Region 

The cusp region as shown in figure 1 is located on the antisolar side of the magnetosphere 
where the dipolar field geometry is transformed into the field geometry of the magnetotail. 
The cusp region is characterized by a dipolar field, high density of charged particles, and an 
abrupt change in the magnetic field intensity at the outer boundary as shown in figure 1 2 
(ref. 14). The region is part of the Van Allen radiation belts and is located between about 8 
to 16 R e and within geomagnetic latitudes of about ±25°. The position of the outer 
boundary of the cusp region varies daily. This boundary is characterized by abrupt and large 
changes in particle density, and has a variation in field intensity of the order of 5 to 20 7 (5 
to 20 nT) (ref. 14). The magnitude of the magnetic field changes from about 30 7 (30 nT) 
on the cusp side of the outer boundary to about 15 7(15 nT) on the magnetotail side of 
this boundary. 


♦Computational model of the geomagnetic field is discussed in appendix C. 
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Figure 9. Magnetic field at synchronous altitude from internal source (Computed from a 
24-term harmonic expansion for 3 November 1968). 
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Figure 10.— ATS-1 magnetic field measurements at synchronous altitude showing diurnal 
variation on 4 January 1967, at 150° west longitude (ref. 13). 
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Figure 11 .— ATS - 1 magnetic field measurements at synchronous altitude showing horizontal 
intensity variations during a magnetic storm (ref. 13 ). 
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Figure 12.— Observed magnetic field intensity and particle density in cusp region (ref. 14). 
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2.1.8 Magnetopause 

The magnetopause as shown in figure 1 is the boundary layer that separates the relatively 
strong geomagnetic field from the magnetosheath. This boundary layer has hourly and daily 
displacements of 1 to 3 R ( but generally extends from about 10 R e at the subsolar side to at 
least 80 R e on the antisolar side of the magnetosphere (ref. 15). During a magnetic storm, 
the changes in the interplanetary magnetic field and in solar wind intensity cause the 
magnetopause to expand or contract at rates greater than 50 km/sec (ref. 13). 

The magnetopause is a barrier region some tens of km in thickness in which the magnetic 
field intensity is low and the magnetic field direction shifts by as much as 180 degrees. The 
magnetopause is the transition zone from field intensities greater than 40 y (40 nT) in the 
magnetosphere to field intensities of about 10 to 30 7 ( 10 to 30 nT) in the magnetosheath 
as shown in figure 13 (ret. 16 ). 



Figure 13. Magnetic field intensity in magnetosheath from Mariner V observations (ref. 16). 
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2.1.9 Magnetosheath 

The magnetosheath, as shown in figure 1, is the transition zone in which the interplanetary 
field and the solar wind are deflected after passing through the shock front. In this region, 
the weak interplanetary field is compressed and deflected. The field compression increases 
the magnetic field intensity in the magnetosheath to a nominal 10 to 30 7 ( 10 to 30 nT) as 
shown in figure 13 (refs. 16 and 17), however, during intense magnetic storms the tield has 
exceeded 100 7 (100 nT). The magnetosheath field intensity is also characterized by rapid 
oscillations with periods of 4 to 15 seconds and with 2 to 57 (2 to 5 nT) peak-to-peak 
amplitudes (ref. 18). 

The thickness of the magnetosheath in the Earth-Sun direction is a nominal 4 R e ; it can 
shrink to about 2 ± 0.5 R e during a magnetic storm that occurs during a year of low solar 
activity. As any magnetic storm activity decreases, the magnetosheath increases to its 
nominal thickness (ref. 15 ). 


2.1.10 Shock Front 

The shock front, as shown in figure 1 . is the outermost boundary of the magnetosheath and 
is created by the interaction of the solar wind with the geomagnetic field. Observations 
made at the shock front indicate the presence of large amplitude magnetic field oscillations 
of 10 to 50 7 (10 to 50 nT) as shown in figure 14 (ref. 18). On the magnetosheath side of 
the shock front, erratic and periodic oscillations occur with amplitudes of 2 to 5 7 (2 to 5 
nT) and periods of about 4 to 15 seconds. At the shock front, higher frequency oscillations 
occur with amplitudes ranging from 5 to 10 7 (5 to 10 nT) and frequencies on the order of 
1 to 3 Hz and higher (ref. 19). On the interplanetary side of the shock front, oscillations 
have been known to occur with amplitudes of 2 to 5 7 (2 to 5 nT) and periods ot about 20 
to 60 seconds (ref. 1 8 ). 


2.1.11 Magnetotail 

The magnetotail, as shown in figure 1, is formed by the solar wind sweeping back the 
geomagnetic field lines in the antisolar direction. Magnetotail field lines from southern 
hemisphere are predominantly directed in an antisolar direction and separated by the 
neutral sheet from field lines from northern hemisphere, which are predominantly solar 
directed. Pioneer VII observations indicate that the magnetotail extends to at least 1000 R e , 
and IMP-1 observations indicate that it has a radius of about 20 R e at a geocentric distance 
of 30 R e (refs. 14 and 20). Direction of the field in the magnetotail is generally close to the 
plane of the ecliptic and at a solar coordinate angle <p gE between 270° and 360° (ref. 2 1 ).* 


♦Solar ecliptic coordinates are defined in appendix A. 
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Figure 14. Observed magnetic field oscillations at shock front (ref. 18). 

The configuration of the magnetotail is well defined to a geocentric distance of about 
80 R e . At distances of 500 to 1000 R e the configuration appears to be characterized by 
separate tubes of magnetic filaments (ref. 22 ). 

The magnetic field intensity in the magnetotail to a geocentric distance of about 8 to 10 R e 
decreases at a rate expected for the Earth’s dipolar field? Beyond 10 R e , Explorer 33 
observations in the magnetotail indicate that the field intensity generally decreases as 
distance from the Earth increases. This decrease is from about 20 7 (20 nT) at 10 R e to 
about 7 7 (7 nT) at 80 R e . Between 1 0 R e and 30 R e , IMP 1 observations indicate that the 
median field intensity is about 16 7 (16 nT) (ref. 22). Ness (ref. 4) indicates that from 
20 R e to 60 R e , a decrease from about 16 to 87 (16 to 8 nT) is to be expected. 
Pioneer VII observations indicate that the field intensity is generally between 4 and 16 7 (4 
and 16 nT) at about 1000 R e . The currents within the neutral sheet cause magnetic Held 
depressions that extend northward and southward from the neutral sheet to a distance of 
about 5 R e and may extend lengthwise along the neutral sheet for distances greater than 

1 5 R e 


*See presentation of dipolar model in appendix C. 
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2.1.12 Neutral Sheet 


The neutral sheet, as shown in figure 1, is a region within the magnetotail where abrupt 
reversal of the magnetic field from the solar to the antisolar direction occurs and in which 
magnetic intensity decreases to a low value. The neutral sheet extends from about 
10±3R e to 1000 R e or to the limit of the magnetotail (ref. 23). The sheet is 
approximately parallel to the solar-magnetospheric equatorial plane* with a position above 
the plane during summer months and below the plane during winter months (ref. 23). 
During a magnetic storm, the front boundary of the neutral sheet has been observed to be 
displaced to geocentric distances beyond 13 R e , returning to about 10 R e as storm activity 
decreases (ref. 24). 

The magnetic field intensity in the neutral sheet is about 1 to 4 7 (1 to 4 nT) on the dawn 
meridian side and less than 1 7 (1 nT) on the side near the noon-midnight meridian plane. 
The thickness of the neutral sheet varies from about 500 km to about 5000 km. 


2.2 Models of the Earth’s Main Field 

The analytical models of the Earth’s main (geomagnetic) field discussed in this monograph 
represent an undistorted, steady-state field and do not account for the effects of external 
sources or localized Earth anomalies. An analytical model that considers some external 
sources has been developed by Mead (ref. 12), but it is a simplified model. 


2.2.1 Spherical Harmonic Expansion Model 

The spherical harmonic expansion model of the geomagnetic field is based on a series 
expansion in spherical coordinates of the geomagnetic potential devised by Gauss in the 
1830’s and which satisfies LaPlace’s and Maxwell’s equations.** The accuracy of this model 
derives from the use of a large number of coefficients in the harmonic expansion. Cain 
Hendricks, Langel, and Hudson (ref. 25) proposed an expansion of 120 coefficients with 
associated first and second time derivatives for use with the model for the “International 
Geomagnetic Reference Field - 1965.” With these coefficients, the model has expected 
maximum errors of the order of a few tens of gammas in total field intensity at the Earth’s 
surface (except for Antarctic regions). At satellite altitudes of about 4 R e , errors are reduced 
by a factor of 3 to 4, but at higher altitudes the effects of external sources increase so that 
the model gradually becomes less representative of the actual environment, particularly 
beyond 6.6 R e (synchronous altitude). 

The spherical harmonic expansion model of Cain, et ai. is currently considered 
the best analytical representation of the geomagnetic field and has been adopted 


* Solar magnetospheric coordinates are defined in appendix A. 

** Analytical presentation of field models is in appendix C. 
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for this monograph. The abbreviated versions of this model discussed in 2.2.2 and 2.2.3 are 
less accurate but often useful for quick estimates of the field. 


2.2.2 Quadrupole Model 

The quadrupole model of the geomagnetic field is derived by using the first eight terms of 
the spherical harmonic expansion model. This model includes both dipolar and quadrupolar 
effects. At geocentric distances greater than 1.2 R e , the maximum error between computed 
and actual fields is about 20% in magnitude and 10° in direction. At distances less than 

1.2 R e , the error increases until it reaches about 40'% in magnitude and 20° in direction at 
the Earth’s surface.* This large error is caused by deviation of the Earth’s surface field from 
a dipole or quadrupole field. Figure 15 (ref. 26) illustrates the magnitude and angular 
deviations of the quadrupole model as a function of geocentric distance. 


2.2.3 Dipole Models 

The spin-axis dipole model of the geomagnetic field is derived by using the first term of the 
spherical harmonic expansion model. This model corresponds to a magnetic dipole located 
at the Earth’s center with its magnetic poles aligned to the Earth’s spin axis. The spin-axis 
dipole model yields field magnitudes that are generally lower than the actual field; this 
model is the least accurate of the models of the geomagnetic field.* 

The tilted dipole model is obtained by using the first three terms of the spherical harmonic 
expansion model. This model has the dipole’s axis tilted away from the Earth’s spin-axis by 
11.44° so that magnetic North Pole is near Thule, Greenland. At geocentric distances greater 
than 3 R e , solar wind and solar plasma effects cause the actual magnetic field to be 
significantly different from the computed tilted dipole field. At 3 R e , maximum errors are 
about 9% in magnitude and 7° in direction. At distances less than 3 R e , magnitude errors 
may be as great as 40 to 60%.* Figure 1 5 illustrates the magnitude and angular deviations of 
the tilted dipole model as a function of geocentric distance. 


2.3 Interplanetary Field 

The interplanetary field is the Sun’s magnetic field which is represented as a sectionalized, 
spiral field. This field extends radially from the solar photosphere for a few Sun radii and 
then begins to form into a spiral as shown in figure 16 (refs. 27, 28 and 29). Beyond several 
solar radii (outside the region indicated by the dashed circle in figure 16) the rotational 
interaction causes the interplanetary field to form into an Archimedes spiral. The field 
polarity varies with the sectors, a sector of one polarity often being as large as 180°. IMP-1 
detected four sectors of alternating polarity as shown in figure 16. The interplanetary field 


* Analytical presentation and accuracy comparison of field models are in appendix C. 
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GEOCENTRIC DISTANCE (R e ) 

Figure 15. Maximum deviations of quadrupole and tilted-dipole models 
from spherical harmonic expansion model (ref. 26). 
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Dashed Circle Indicates 
Beginning of Spiral Effect 
at Several Sun Radii. 


^Figure 16. Sample structure of interplanetary field (refs. 27, 28, and 29). 


intersects the Earth’s orbit in the plane of the solar equator (about 7° from plane of 
ecliptic) at an angle of about 45° to the Earth-Sun line as illustrated in figure 16 when the 
solar wind velocity is 450 km/sec (ref. 28). Microstructure in the interplanetary field is the 
source of magnetic fluctuations of less than one-hour duration which are characterized by 
shock waves and plasma-magnetic field discontinuities. Mesostructure field effects produce 
magnetic variations of about 1 to 100 hours and are characterized by magnetic filaments 
(tubular groupings of lines of force) and low-frequency magnetohydrodynamic waves. 
Macrostructure field effects produce magnetic variations of over 100-hour duration and are 
characterized by sectoring and very wide magnetic filaments (ref. 30). as illustrated in figure 
16, with diameters of an order of 3X10 6 km. Pioneer VI observations provide evidence of 
the existence of magnetic filaments and also indicate that they are sometimes intertwined 
and corotate with Sun. The rotating sector structure and polarity of the interplanetary field 
has been increasingly studied in recent years. Satellite measurements in the vicinity of 1 
A.U. have been extrapolated in the ecliptic plane to obtain mapping of the field between 0.4 
and 1.2 A.U. (ref. 3 1 ). A solar flare can send an intense magnetic disturbance propagating 
outward into the interplanetary field disrupting its structure. After the disturbance, a new 
structure forms that is very similar to the preceding one. 

The magnitude of the interplanetary field in the vicinity of the Earth is fairly constant with 
a nominal value of about 6 y (6 nT) and a range of 2 to 40 7 (2 to 40 nT) (ref. 30). At the 
Sun’s equator and outside of active centers the solar magnetic field intensity is estimated at 
about 40 gauss (.0004 T). An average interplanetary intensity of 2 7 (2 nT) implies that the 
magnetically active region of the photosphere has a magnetic field intensity of about 100 
gauss (.01 T) (ref. 29). 

Large amplitude disturbances of the interplanetary field are associated with transient solar 
features which occur in centers of local activity. Strong magnetic fields are characteristic of 
sunspots, occasionally reaching 3000 gauss in a center (ref. 15). The polarity of the lead 
spots (in the sense of rotation, that is, west) reverse with an 1 1-year period and the sun’s 
main field reverses with 22-year periodicity (ref. 32). 


2.4 Magnetic Fields of Planets and Their Natural Satellites 

Spacecraft magnetometer observations have provided upper bounds on the dipole moments 
of Venus, Mars, and the Moon. Radio observations show that Jupiter has a strong magnetic 
field, and modulation of decameter emissions indicate that its satellite Io may also have 
magnetic effects (ref. 33). There have been no observations yielding positive evidence for 
existence of magnetic fields for the planets Mercury, Saturn, Uranus, Neptune, and Pluto 
nor any associated natural satellites. 


2.4.1 Venus 

Spacecraft observations indicate either that Venus has no magnetic field or that it has a field 
of the same order as the interplanetary field. Mariner V observations indicate an 
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interplanetary field intensity of about 8 7 (8 nT) near the orbit of Venus and an upper limit of 
the magnetic dipole moment for Venus of about lO 3 times the Earth’s dipole moment*, 
implying that Venus’ magnetic field at the equator is less than 35 7 (35 nT). Figure 17 
(ref. 16) illustrates the interaction of the solar wind and Venus. In this concept, the solar 
wind plasma flows around Venus leaving a cavity behind the planet. Mariner V observations 
indicate a shock front bounding a region similar to Earth’s magnetosheath. The solar wind 
plasma behind the shock front is bounded by an anemopause which is about 50 km above 
the surface. The anemopause is the interaction layer between Venus’ atmosphere and the 
solar wind. The expansion fan illustrated in figure 17 could be formed by solar plasma 
expanding into the cavity behind Venus with a rarefaction expansion wave moving away 
from the cavity in reaction to the inward plasma flow (refs. 16 and 34). 


SOLAR WIND 


SOLAR WIND 



Figure 17.— Interaction model of solar wind and Venus (ref. 16). 
*7.99 X 10 3 * pole -cm (10.04 X 10*^ weber-meter) is Earth’s magnetic moment. 
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2.4.2 Mars 


Observations made by Mariner IV at 3.9 Mars radii from the center of the planet indicate 
that no definite effect could be attributed to a Martian magnetic field. Citing these 
observations, Smith (ref. 35) estimates the Martian magnetic moment to be less than 10 4 
times Earth’s magnetic moment.* A magnetic moment of that order implies a Mars surface 
magnetic field of about lOOy (100 nT). If Mars has a negligible magnetic field, the solar 
wind interaction with Mars may be similar to the Venus interaction. However, the magnetic 
field of Mars may be strong enough to form a magnetosphere with a structure similar to the 
Earth’s magnetosphere (ref. 16). 


2.4.3 Jupiter 

Jupiter is expected to have a magnetosphere very similar to that of the Earth and extending 
about 40 to 50 Jupiter radii along Sun-Jupiter line as shown'in figure 18 (refs. 16 and 36). 
Observations of radio emissions from the planet indicate that its magnetic field is at least 10 
times as intense as Earth’s (ref 4). Smith (ref. 16) reports the estimated value of the surface 
field to be about 5 gauss (.0005 T). Jupiter’s magnetic field resembles a magnetic dipole 
field with its axis apparently displaced off the planetary center and tilted at 9° or 10° to 
the axis of rotation; the magnetic polarity is opposite that of Earth’s dipole field (ref. 37). 
Because of a short rotation period (less than 10 hours) conjunctive with a suspected massive, 
fluid, conducting center, Jupiter’s magnetic moment is expected to be greater than the 
Earth’s by as much as a factor of 1 0 4 * (ref. 16). 


2.4.4 Saturn 

Observations of possible decameter radio emissions from Saturn indicate that this planet 
may have a moderately strong magnetic dipole moment (ref. 38). In every case the weakness 
and uncertainty of the radiation defy positive identification. Thus if Saturn is a source of 
nonthermal radio emissions, the output is far less frequent and much less intense than that 
from Jupiter. Further observations are required, but Saturn, with its similarity in size and 
appearance to Jupiter, may have a strong magnetic field whose configuration may be 
influenced by the presence of Saturn’s rings (ref. 39). 


2.4.5 Moon 

Measurements indicate that Moon’s surface magnetic field is less than 8 7 (8 nT) and 
magnetic moment is less than 10 s times the Earth’s magnetic moment* (ref. 40). Explorer 
35 data indicates that a lunar shock front does not exist and that the interplanetary field is 
not significantly distorted by the Moon as shown in figure 19 (ref. 16). Without a shock 
front, the Moon probably absorbs the solar wind plasma that strikes it with very few 
particles being deflected. 


‘7.99 x io 2s 


pole-cm (10.04 X 10 16 


weber-meter) is Earth’s magnetic moment. 
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INTERPLANETARY 
FIELD - 0.7r 



Figure 18, Magnetosphere of Jupiter (refs. 16 and 33); an idealized 
schematic cross-section in the ecliptic plane. 

The Moon’s wake region is characterized by a solar wind umbra and penumbra. Wang 
(ref. 41) predicts a long non-axial, symmetrical wake behind the Moon (the plane of 
symmetry is defined by the local solar wind velocity and the interplanetary field direction). 
The interplanetary field in the solar wind umbra increases about 3 7 (3 nT) with a relative 
decrease of about 4 7 (4 nT) observed on either side of the increase. The penumbra is 
characterized by decreases in the interplanetary field which are externally bounded by 
increases in field intensity (ref. 41). 

RAREFACTION 



INTERPLANETARY FIELD 


Figure 19. Interaction model of solar wind and Moon (ref. 16). 
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3. CRITERIA 


The analytical models and descriptions of the magnetic Helds presented here should be used 
in the design of space vehicles, space vehicle systems, and the experiments and 
instrumentation for space applications. 


3.1 Earth’s Surface 

The surface geomagnetic field for the geographical region of interest may be obtained from 
magnetic charts issued by the U.S. Naval Oceanographic Office, Washington, D.C. 20390 
(ref. 8). Figure 3 is a simplified version ot one such chart showing the total magnetic field 
intensity at the Earth’s surface. 


3.2 Earth’s Surface to Synchronous Altitude 

The magnetic field environment from the Earth’s surface to synchronous altitude 
(geocentric distance of 6.6 R e ) can be described for most design purposes by the spherical 
harmonic expansion model (appendix C). 


3.2.1 Secular Variations 

Secular variations can be as large as 150 7 (150 nT) per year. These variations are included 
in the spherical harmonic expansion model (appendix C). 


3.2.2 Temporal Variations 

Temporal variations are not included in the spherical harmonic expansion model but should 
be considered if these variations could significantly affect a spacecraft’s instrumentation or 
scientific experiments. Expected field changes for temporal variations are as follows. Diurnal 
variations are generally less than 200 7 (200 nT). Variations caused by magnetic storms 
vary from about a 50 7 (50 nT) increase at the surface during the sudden 
commencement phase to a greater than 400r (400 nT) decrease at the surface during main 
phase (the time sequences and durations of the phases are described in 2.1.3). Sudden 
impulse variations are ±18 7 (±18 nT) maximum at the surface. Micropulsations cause a 
variation in the field of a few 7 (nT) at the surface. 


3.2.3 Ionospheric Currents 

Magnetic field changes caused by ionospheric currents are less than 2000 7 (2000 nT) and 

do not make significant contributions to the main geomagnetic field except at .geocentric 
distances below about 2 to 3 R e . 
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3,2.3. 1 Auroral Currents 

Auroral currents (polar substorm) can cause field changes as large as 2000 y (2000 nT), but 
the average is about 1 500 y ( 1 500 nT) with a duration of 30 to 1 20 minutes. 


3. 2. 3. 2 Solar Currents 

At the Earth’s surface, S q currents can cause field changes in the middle latitudes of 20 to 
40 7 (20 to 40 nT). The magnitude of this change increases to 100 to 200 y (100 to 200 nT) 
near the magnetic equator. 


3. 2. 3. 3 Lunar Currents 

At the Earth’s surface, lunar currents cause insignificant changes of about 1 7 (1 nT) with a 
12 hour period. 


3.3 Synchronous Altitude 

At synchronous altitude, an average magnetic field of about 138 7 (138 nT) is to be 
expected. Approximations to the magnetic field in this region may be obtained from the 
spherical harmonic expansion model (appendix C), but expansion beyond the first 33 terms 
in this model is not warranted because contributions from external sources are not included. 
The field at synchronous altitude is highly variable during a magnetic disturbance with 
increases of 50 percent, decreases to zero, and even reversals of direction (see 2.1.6). 


3.4 Beyond Synchronous Altitude 

At distances between 6.6 and about 10 R e , the spherical harmonic expansion model 
(appendix C) will give an approximation to average quiet day conditions but diurnal and 
temporal effects in this region are of major significance. Because of the limited applicability 
of the spherical harmonic model in this region, computations using more than about 20 
terms are not warranted. About 10 R e and beyond, distortions in the geomagnetic field 
cause deviations from the spherical harmonic expansion model; therefore, estimates of the 
magnetic field are given for the following regions. 


3.4.1 Cusp 

The cusp region is part of the Van Allen radiation belts and is located on the night side of 
Earth between about 8 to 16 R e and within geomagnetic latitudes of about ±25°. The cusp 
region is characterized by dipolar field geometry and an abrupt and large change in particle 
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density and field intensity at its outer boundary. The magnitude of the magnetie field 
throughout the cusp region is not well known, but it changes from about 30 7 (30 nT) near 
its outer boundary to about 15 7(15 nT) on the magnetotail side of this boundary. 


3.4.2 Magnetopause 

The magnetopause is the transition zone from field intensities greater than 40 7 (40 nT) in 
the magnetosphere to field intensities of about 10 to 30 7 (10 to 30 nT) in the 
magnetosheath. This transition zone or boundary layer is about 100 km thick and has 
hourly and daily displacements of as much as 1 to 3 R e but generally bounds the 
magnetosphere from about 1 0 Ron the subsolar side to at least 80 R on the antisolar side. 

“ C 


3.4.3 Magnetosheath 

The magnetic field intensity in the magnetosheath is a nominal 10 to 30 7 (10 to 30 nT), 
but during intense magnetic storms, the field intensity has exceeded 100 7 (100 nT). The 
magnetosheath is characterized by rapid field oscillations with periods of 4 to 1 5 seconds 
and with 2 to 5 7 (2 to 5 nT) peak-to-peak amplitudes as well as by periodic and incoherent 
noise components. The thickness of the magnetosheath in the Earth-Sun direction is a 
nominal 4 R e but can shrink to about 2 ± 0.5 R e during a magnetic storm that occurs during 
a year of low solar activity. 


3.4.4 Shock Front 

Observations at the shock front indicate the presence of large amplitude magnetic field 
oscillations of 10 to 50 7 (10 to 50 nT). On the magnetosheath side of the shock front, 
erratic and periodic oscillations occur with amplitudes of 2 to 5 7 (2 to 5 nT) and periods 
of about 4 to 1 5 seconds. At the shock front, higher frequency oscillations occur with 
amplitudes ranging from 5 to 10 7 (5 to 10 nT) and frequencies on the order of 1 to 3 Hz 
and higher. On the interplanetary side of the shock front, oscillations occur with amplitudes 
of 2 to 5 7 (2 to 5 nT) and periods of 20 to 60 seconds. 


3.4.5 Magnetotail 

The magnetic field intensity in the magnetotail to a geocentric distance of about 8 to 10 R e 
decreases at a rate expected for the Earth’s dipolar field. From about 10 to 80 R e , the field 
decreases from about 20 7 (20 nT) to about 7 7 (7 nT). Between 10 R e and 30 R e , the 
median field intensity is about 16 7 (16 nT). From 20 R e to 60 R e , a decrease from about 
16 to 8 7 (16 to 8 nT) is to be expected. At about 1000 R e , the field intensity is generally 
between 4 and 16 7 (4 and 16 nT). The direction of the field in the magnetotail is generally 
close to the plane of the ecliptic and at a solar coordinate angle <f> SE (appendix A) between 
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270° and 3b0°. Held lines Irom the Earth's southern hemisphere are predominantly 
directed in an antisolar direction and separated by the neutral sheet from field lines from 
the northern hemisphere which are predominantly solar directed. 


3.4.6 Neutral Sheet 

The neutral sheet is a region within the magnetotail where abrupt reversal of the magnetic 
field from the solar to the antisolar direction occurs and in which magnetic intensity 
decreases to 4 7 (4 nT) or less. The magnetic field intensity in the neutral sheet is about I 
to 4 7 (1 to 4 nT) on the dawn meridian side and less than I 7 (1 nT) on the side near the 
noon-midnight meridian plane. This sheet extends from about 10 ± 3 R e to 1000 R e or to 
the limit of the magnetotail. Its thickness varies from about 500 km to about 5000 km. 


3.5 Interplanetary 

An analytical model is not available for use in evaluating the interplanetary magnetic field. 
This field can be described as having a spiral type geometry divided into sections of constant 
polarity. The interplanetary field has a range of 2 to 40 7 (2 to 40 nT) with an average 
magnitude of 67(6 11 T) near Earth s orbit. 


3.6 Planets and Their Natural Satellites 

Analytical models are not available for any of the planets (other than Earth) or their natural 
satellites Space flight data has established bounds for the magnetic fields of Venus as less 
than 35 7 (35 nT), Mars as less than 100 7 ( 100 nT) at the surface, and the Moon as less 
than 8 7 (8 nT) at the surface. The surface magnetic field at the equator of Jupiter is 
estimated to be about 5 gauss (.0005T). Jupiter’s dipole moment is expected to be 10 
times that of the Earth’s (ref. 16). In the vicinity of Jupiter’s shock front, the interplanetary 
field is estimated to be about 0.7 7 (0.7 nT). 
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APPENDIX A 

COORDINATE SYSTEMS AND NOMENCLATURE 


Several coordinate systems are commonly used to define magnetic quantities in the context 
of convenient geometry for the geomagnetic field, the magnetosphere, and the solar system. 
The following coordinate systems are defined in this appendix: geocentric, geodetic, 
geomagnetic, local geomagnetic, solar ecliptic, and solar magnetospheric. Because geographic 
and geodetic coordinate systems are equivalent, the terms geographic and geodetic are used 
synonymously. 
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Geocentric Coordinate System 


z E 



Figure A-l . Geocentric and geodetic coordinate systems. 


SYMBOLS 

X E Positively directed along intersection of Greenwich 

Meridian and Equatorial planes 

Z E Positively directed north along Earth’s mean spin axis 

Y e Z e xX e 

X Geocentric longitude, measured positive eastward from 

Greenwich meridian to local meridian 
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5 Geocentric latitude, declination measured positive north 

(geocentric and geodetic latitude differs by 1 1.6 minutes 
maximum at 45° latitude.) 

0 Geocentric colatitude = 90° - 5 

R e Mean Earth radius, 6371.2 km 

X Geodetic longitude = geocentric longitude 

<p Geodetic latitude, measured positive north from 

Equatorial plane to normal to Spheroid 
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Geomagnetic Dipole Coordinate System 


240 ° 270 ° 300 ° 330 ° 0 ° 30 ° 60 ° 90 ° 120 ° 150 ° 180 ° 210 ° 240 ° 



Figure A-2 - Geomagnetic dipole coordinate system (ref. 9). 
SYMBOLS 

A Dipole eastward longitude, angular distance from 

great circle of spherical Earth which passes through 
both geographical and calculated dipole poles. 

<I> Dipole latitude, angular distance north or south of 

dipole equator. 
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DEFINITIONS 


Dipole North Pole. Defined by axis of calculated central dipole whose magnetic field is 
best fit to main geomagnetic field over Earth’s surface and whose axis passes through 
geocenter. Defined to be at 78.5°N, 69.0°W in northwestern Greenland. 

Dipole Equator. -Defined as great circle of spherical Earth which is normal to central dipole 
axis. 

Dip North Pole. -Shifting point on Earth’s surface where geomagnetic field lines are vertical. 
Location changes by a few kilometers during any day in response to transient magnetic 
fields. Located northwest of Hudson Bay at about 73° N, 98°W. 

Geodetic to Dipole Transformation - This is defined by the following equations (ref. 11): 


. cos (p sin (X - 291°) 

sin A = 

cos <1> 


sin < f> = sin <p cos 11.7° + cos 0 sin 1 1.7° cos (X - 291°) 
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Local Geomagnetic Coordinate System 



Figure A4. Local geomagnetic coordinate system (ref. 9). 

SYMBOLS 

B Total computed magnetic field intensity (total magnetic 

flux density - SI nomenclature) 

F Total measured field intensity 

X(B x ) North component of B 
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Z (B Z ,V) 
Y (B y ) 


Vertical component of B 


East component of B (Z x X) 


1 Inclination or dip angle 

H Horizontal field intensity (B cos I) 

D Declination angle 


Solar Ecliptic Coordinate System 



Figure A-5.— Solar ecliptic coordinate system (ref. 42). 


SYMBOLS 


^SE 

Positively directed 

^SE 

Positively directed 

^SE 

Z x x 

E S E 


toward Sun along Earth-Sun line 
toward north ecliptic pole 
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Solar ecliptic latitude 


^SE 

<p SE Solar ecliptic longitude 

Solar Magnetospheric Coordinate System 



Figure A-6.— Solar magnetosphiric coordinate system (ref. 42). 

SYMBOLS 

M e Positively directed northward parallel to Earth’s 

magnetic dipole axis 

S E Positively directed northward parallel to Earth’s 

spin axis 

X SM Positively directed toward Sun along Earth-Sun 

line (parallel to X SE ) 

^SM ^SM* ^SM 

Y sm Parallel to M f; x X sm 
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APPENDIX B 

DIMENSIONS, UNITS, AND CONVERSION FACTORS 


Several systems of units have been commonly used to describe the characteristics of 
magnetic fields. The system of electrostatic units (esu) uses Coulomb’s law; and the system 
of electromagnetic units (emu) uses the law of attraction between currents. The Gaussian 
system expresses magnetic quantities in emu and electric quantities in esu. In the Gaussian 
system, the magnetic flux density (B) and the magnetic field intensity (H) are used 
interchangeably. This interchangeability is based on the fact that B = juH. where the 
permeability ju, is unity for a non-magnetic medium, in which case B = H. The unit for B and 
H was the gauss until 1932. After 1932, H was expressed in oersted. Another unit used for B 
and H is gamma ( 7 ) which was first introduced in studies of geomagnetism. When fi is unity, 
gamma is equal to ICC 5 gauss and is validly used for both field intensity and flux density 
(ref. 43). In 1960, the Eleventh General Conference on Weights and Measures adopted the 
metric system of units (meter, kilogram, second, ampere. Kelvin degree, and candela) as the 
international standard. This system was named the International System of Units (SI). In SI 
units, the tesla ( I ) is the unit of magnetic flux density and the ampere per meter ( A/m) is 
the unit of magnetic field strength (ref. 44). 

Dimensions and units used in this monograph are in terms of both the Gaussian system and 
the International System. Relations between the SI metric system of units and emu and esu 
systems are given below for the following quantities: magnetic field intensity, magnetic flux, 
magnetic field strength, magnetic dipole moment, and current (ref. 44). 
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TABLE B - 1 Magnetic Field Intensity (Magnetic 
Flux Density in SI Nomenclature)* 


System 

of 

Units 

equals 

one 

T 

r 

Y 

SI 

webers / meter 2 ( Wb / m 2 ) 
or tesla (T) 

1 

10 4 

10 g 

emu 

gauss ( r ) 

10' 4 

1 

10 s 

** 

gamma ( Y ) 

10-9*** 

M 

1 


* Vector components of field intensity are defined 
in figure A -4 of appendix A. 

** Not part of a system. 

*** 1 nanotesla (nT) 


TABLE B-2 Magnetic Flux (4>^) 


System 

of 

Units 

N. equals 

one 

Wb 

line 

SI 

weber (Wb) 

1 

10* 

emu 

maxwell (line) 

10' 8 

1 


TABLE B-3 Magnetic Field Strength (H) 


System 

of 

Units 

equals 

one 

A/m 

Oe 

Y 

SI 

ampere / meter ( A / m ) 

1 

4 it x 10' 3 

4 rrX 10 2 

i 

emu 

oersted ( Oe ) 

10 3 / 4fr 

1 

10 s 

# 

gamma (y) 

10' 2 I4 n 

mm 

1 


* Not part of a system 
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TABLE B-4 Magnetic Dipole Moment (m) 


System 

of 

Units 

equals 

one 

1 

1 


Loop: 

polecm 

Dipole: 

polecm 

SI 

amperemeter 2 (A-m 2 ) 

1 

4 rr xIO 7 

10 3 

10 3 

weber-meter (Wb*m) 

10 7 

1 

io 10 

IO 10 

4 TT 

4 TT 

4tt 

emu 

pole - centimeter ( pole - cm 
or upc ) 

10 3 

4tt x10“’° 

1 

1 


Note: In cgs units the following are often used interchangeably for magnetic 
moment: pole -cm; gauss -cm 3 ; dyne - cm / oersted; dyne - cm / gauss. 


TABLE B -5 Current 


System 

of 

Units 

equals 

one 

A 

abA 

statA 

SI 

ampere ( A ) 

1 

0.1 

3x10 9 

emu 

abampere ( abA ) 

10 

1 

3X10'° 

esu 

statampere ( statA ) 

- xIO 9 
3 

I X 10 10 

3 

1 
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APPENDIX C 

MODELS OF THE EARTH’S MAGNETIC ENVIRONMENT 


Spherical Harmonic Expansion Model 


The spherical harmonic expansion model is a series expansion in spherical coordinates of the 
magnetic flux density (B) which is defined as the negative gradient of the geomagnetic 
potential function V, where 


oo n 

"-"•EE 

n = 1 m - 0 


K\ n+l 
R / 


(g™ cos mX + h™ sin mX (cos d) (C-l) 


R e = mean Earth radius (6371.2 km) 

R = distance from geocenter 
g m and h m = harmonic coefficients 

n n 

X = geocentric east longitude 

6 = geocentric colatitude (90° - 6) 

5 = geocentric latitude 

P™ (cos 6) = Schmidt function - partially normalized spherical 

harmonic function of degree n and order m with 
m^n. 

The harmonic coefficients g n and h™ are determined by analyzing measurements of the 
geomagnetic field components made at a large number of points on the Earth’s surface. 
Table C-l (ref. 25) gives a set of 120 of these coefficients with their first and second 
derivatives. 
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The Schmidt functions are defined as: 


P n , m (cos0) 


for m = 0 


P n (cos 6 ) =. 


2 ( n - m ) ! 


( n + m ) ! 


P„,« (cosO) 


for m 4- 0 


(C-2) 


The associated Legendre function, P n m (cos 9). is given by 


P„ m (cos 0) = 


(2n) ! 


sin m 0 


2 n n ! (n - m) ! 


(cos 0) n - m 


(n - m) (n-m - 1) 
2 (2n- 1) 


( cos 9 ) 


n -m - 2 


(n-m) (n-m-1) (n-m-2) (n-m-3) 
2(4) (2n- 1 ) (2n - 3 ) 


(COS 0)"-m-4 + . . . 


(C-3) 


To calculate components of B at some point in space, the following relations are used: 


oo n m 

1 3V V^V^/ R e\ n+2 dP n < COS °) , n. , . .m ■ ,, , ri , 

g v = = S > I J 2, [g cos mA + h sin mA| (C-4) 

n = 1 m 0 R 


R 3 0 


d0 


oo OO 

1 3 V R e\ mP (cost?) m . ,. , rc . 

= - > y [ *] l l (g sin mX - h cos m\] (C-5) 


Rsin 6 3X ^ 

n = 1 m - 0 


sin 0 


R z 


3 V 
3 R 


00 . n / R \n + 2 

^ / — Ij (n+l)P™ (cos 6 ) | g' n cos mA + h m sin mA] (C-6) 


n = 1 m - 0 


♦Refer to appendix A for definition of geocentric coordinate system. 
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where: 


dP n (cost?) 

d e 


— (n+l) 

2 


'/2 


P (cos 0 ) 


for in = 0 


and 


dP n (cos (?) 

dF 


- l/. 


'm- 1 


(n + m) (n - m + 1 ) 


'/? 


p' n ' (cos fl) 


(C-7) 


1 m t l . „ 

- '/> (n + H1+1 ) (n-m) P (cos0) tor in =£ 0 


i = 2 

5m- 1 = 1 


for in = 1 
for m 2r 2 


The coordinates of a point or points in space and the number of coefficients (corrected for 
secular variation) needed to obtain the desired accuracy are substituted into equations C-4 
through C-6. Secular variation correction is computed by adding the secular change in the 
coefficient over a period of time to the original coefficient. The corrected coefficient is 
defined as: 


g m (t) = g m (t) + (t-t) g m 

n n n 

h m (t) = h m (t) + (t-t) h m 

n n n 


(t) 

(C-8) 

(t) 

(C-9) 


where t is the epoch in years for which the coefficients were derived and t is the particular 
time of interest (table C-l ). Figure C-l illustrates the relation between the number of 
coefficients (k) used in the equations and the accuracy obtained for k to 63 terms based on 
a set of coefficients published in 1965 (ref. 45). The decimal residual (R) is defined as the 
ratio of the rms of the computed values minus observed values to the observed values: 




(C-10) 
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Table C-l Spherical Harmonic Coefficients (Epoch 1965.0) (ref. 25) 


n 

m 

g 1 V | 

h 

g ( V / yr ) 

h 

g ( V / yr 2 ) 

h 

1 

0 

-30401 2 

000.0 

14.03 

0.00 

- 0.062 

0 000 

1 

1 

- 2163.8 

5778.2 

8.76 

- 3.71 

0.1 14 

- 0.043 

2 

0 

- LS40.1 

OOOO 0 

- 23.29 

0.00 

- 0.154 

0 000 

2 

1 

2997 9 

- 1932.0 

- 0.09 

- 14.31 

-0018 

0.054 

2 

2 

1590.3 

202.9 

- 4.56 

-16 62 

- 0.253 

- 0.016 

3 

0 

1 307. 1 

000.0 

- 8.93 

00 00 

- 0.123 

0 000 

3 

1 

- 1988.9 

-425 4 

- 10.62 

5 20 

-0 027 

0 095 

3 

2 

1 276.8 

227.8 

2.31 

2.53 

0028 

- 0.007 

3 

3 

881.2 

- 133.8 

- 5.89 

- 6.98 

- 0.183 

0.079 

4 

0 

949.3 

000.0 

1.45 

000 

0 00 ! 

0.000 

4 

1 

803.5 

160.3 

0.90 

- 2.19 

-0044 

0 004 

4 

2 

502.9 

- 274.3 

- 1.75 

- 0.14 

0017 

0 056 

4 

3 

- 397.7 

2.3 

0.66 

1.88 

0.007 

-0 035 

4 

4 

266.5 

- 246.6 

-3 01 

-6 52 

-0097 

-0 047 

5 

0 

- 233.5 

000.0 

1.61 

0.00 

0.045 

0.000 

5 

1 

355.7 

5.1 

0.60 

2.24 

0.001 

-0 046 

5 

2 

228.4 

117.8 

3.34 

1.59 

0.075 

0.007 

5 

3 

-28 8 

- 1 14.8 

- 0.04 

- 2.61 

0.008 

- 0.007 

5 

4 

- 157.9 

- 108.9 

- 0.60 

0.50 

0.015 

0 001 

5 

5 

- 62.2 

82.4 

1.76 

- 0.12 

0.056 

- 0.024 

6 

0 

49.2 

00.0 

- 0.42 

0.00 

- 0.006 

0 000 

6 

1 

57.5 

— 12.1 

3.82 

0.05 

0015 

0.020 

6 

2 

- 0.8 

104.4 

0.82 

0.09 

0.0 10 

- 0.01 1 

6 

3 

- 238.3 

56.6 

2.35 

2.55 

0050 

0.015 

6 

4 

- 1.5 

- 23.4 

0.83 

- 1.19 

- 0.011 

-0 029 

6 

5 

- 2.0 

- 14.8 

0.01 

0.33 

0.026 

0029 

6 

6 

- 108.4 

- 13.3 

0.23 

0 84 

0.023 

- 0.010 

7 

0 

72.2 

00.0 

- 0.57 

0.00 

-0014 

0 000 

7 

l 

- 53.7 

- 53.7 

- 0.34 

- 0.96 

- 0.006 

- 0.014 

7 

2 

7.9 

27.4 

- 1.44 

0.01 

-0 034 

0 016 

7 

3 

15.6 

- 8.1 

- 0.90 

0.43 

- 0.004 

0.014 

7 

4 

- 24.3 

7.0 

0.03 

0.75 

- 0.006 

0.005 

7 

5 

- 3.6 

24.3 

- 0.60 

- 0.33 

- 0.027 

- 0.008 

7 

6 

15.5 

- 22.5 

- 0.17 

0.49 

- 0.001 

0.016 

7 

7 

3.6 

-21 4 

- 0.64 

0.90 

- 0.004 

0.011 

8 

0 

8.5 

00.0 

0.35 

0.00 

0.006 

0.000 

8 

1 

6.5 

5.4 

0.50 

- 0.50 

0.008 

- 0.015 

8 

2 

- 9.3 

- 11.7 

1.70 

- 0.21 

0.039 

- 0.012 

8 

3 

- 9.6 

4.2 

- 0.11 

0.03 

- 0.008 

0.005 

8 

4 

- 6.1 

- 15.3 

0.34 

- 0.79 

0.015 

-0 011 

8 

5 

5.5 

4.6 

- 0.07 

0.05 

- 0.002 

- 0.000 

8 

6 

- 8.1 

21.9 

0.43 

0.10 

0.005 

-0003 

8 

7 

130 

- 0.7 

- 0.15 

- 0.36 

- 0.008 

-0 009 

8 

8 

7.4 

-17 1 

- 0.42 

- 0.43 

- 0.007 

-0 003 

9 

0 

10.4 

00.0 

- 0.10 

0.00 

- 0.005 

0000 

9 

1 

5.8 

- 22,4 

- 0.13 

0.66 

- 0.001 

0.022 

9 

2 

7.5 

13 8 

- 1.20 

0.54 

- 0.027 

0.007 

9 

3 

- 15.1 

6.3 

0.08 

0.03 

0.005 

- 0.002 

9 

4 

12.1 

- 3.0 

- 0.08 

0.35 

- 0.007 

0.009 

9 

5 

4.7 

- 1.9 

- 0.39 

- 0.03 

- 0.006 

0.006 

9 

6 

0.2 

9.0 

- 0.36 

- 0.01 

- 0.009 

- 0.001 

9 

7 

1.6 

11.5 

0.47 

0.45 

0.006 

0.009 

9 

8 

0.9 

0.1 

0.37 

“ 0.05 

0.005 

- 0.004 

9 

9 

0.2 

- 1.5 

- 0.48 

0.75 

- 0.009 

0.019 

10 

0 

- 2.9 

0.0 

“ 0.01 

0.00 

- 0.003 

0.000 

10 

1 

" 0.9 

- 0.1 

- 0.13 

- 0.61 

- 0.003 

- 0.012 

10 

2 

- 2.2 

4.5 

0.88 

- 0.64 

0.020 

- 0.014 

10 

3 

0.8 

- 1.0 

- 0.18 

0.02 

- 0.008 

0.001 

10 

4 

- 2.8 

2.6 

0.17 

0.05 

0.007 

0.001 

10 

5 

6.4 

- 4.4 

-0 02 

- 0.63 

0.001 

- 0.011 

10 

6 

4.7 

- 1.3 

0.05 

- 0.07 

0.001 

- 0.001 

10 

7 

- 0.2 

- 3.6 

0.17 

0.07 

0.001 

0.001 

10 

8 

1.8 

4.0 

0.16 

- 0.03 

0.005 

- 0.001 

10 

9 

2.0 

1.0 

0 31 

- 0.02 

0.004 

0.001 

10 

10 

1.1 

- 2.0 

- 0.23 

- 0.45 

- 0.002 

- 0.006 
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where 


I = total number of observations 

F. = scalar field computed at each data point using 
indicated number of coefficients 

F" = respective observed value of scalar field intensity. 

An alternative formulation of this model employs Gauss normalized (Gauss-LaPlace) 
functions, P nm (0), rather than the Schmidt functions, P n m (0). in this case, the magnetic 
flux density is equated to the positive gradient of the geomagnetic potential function 
(B = +VV) (ref. 43 ). 

Computer programs which provide numerical solutions' of the field equations C-4 through 
C-6 are available at NASA Centers as well as most of the major aerospace contractors. 
Generally, these programs allow for variation in the number of coefficients used and thus 
are applicable when limitations in the storage capacity of a particular computer do not 
permit utilization of all the coefficients. 

Reference 43 presents a set of programs which will evaluate the geomagnetic field 
components using any available set of coefficients at any point in normal geodetic 
coordinates (latitude, longitude, altitude). These programs accommodate either Gauss 
normalized or Schmidt functions. 


Quadrupole Model 


The quadrupole model uses the first eight terms of the spherical harmonic expansion to 
obtain an approximation of the geomagnetic field. The geomagnetic potential is 
approximated by: 


V = FT 


(^~~j ( cos O + (gj cos X + hj sin X ) P ( (cost?) 


+ R„ 


R / 


g^ P (cos 0 ) + (g cos X + h sin X) P 2 (cos 0) 


+ (g cos 2 X + h sin 2 X > P (cos 0) 


(C-ll) 
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To calculate the components of B at some point in space, the following are used:* 


1 dV / R e \ 3 o . , 

B x = | — ) [-g sin 0 + cos 0 (g cos X + h sin X) | 

R 3 0 \ R / 1 1 1 


R e \ 4 
R 


-3 sin 0 cos 8 + yj 3 (2 cos 2 9 - 1 ) (g^ cos X + l/ sin X) 


+ (y/~3 sin 0 cos 0) (g 2 cos2X+h 2 sin 2 X) 


(C-12) 


By = 


1 


av _ / R e \ 3 


R sin 6 3X \ R 


g sin X - h cos X 




v/T cos 0 (g l sin X - h 1 cos X ) 


+ sin 0 (g 2 sin 2 X - h 2 cos 2 X ) 
2 2 


(C-13) 


B, 


a v 

a r 


-) 

R 


2 g° cos 0 + 2 sin 0 (g 1 cos X + h 1 sin X) 


M 4 

R / 


— cos 2 0 )(s° ) + s ‘ n ® cos 0 (S* eos X + h 1 sin X) 


2 


2 /V 2 


+ j — sin 2 0j fg 2 cos 2 X + h 2 sin 2 X 


(C-14) 


The accuracy of the quadrupole model compared to the spherical harmonic expansion 
model is illustrated in figures C-2 through C-6 (ref. 26). The magnitude ratio is defined as 
the ratio of the B computed by the approximation model to the B computed from the 
spherical harmonic expansion model. The angular deviation is defined as the field angle 
computed by the approximation model less the field angle computed by the spherical 
harmonic expansion model. 


•Refer to appendix A for definition of geocentric coordinate system. 
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Dipole Models 


Spin-Axis Dipole Model 


The spin-axis dipole model uses the first term of the spherical harmonic expansion to obtain 
an approximation of the geomagnetic field which is generally lower than the actual field. 
The geomagnetic potential is approximated by: 


V = 



g° P° (cosO) 


(C-15) 


To calculate the components of B at some point in space, the following relations are used:* 


B 


X 



0 • n 

g^ sin 0 


(C-l 6) 


B z 


0 


, (A)' 


(C-l 7) 
(C-l 8) 


The accuracy of the spin-axis dipole model compared to the spherical harmonic expansion 
model is illustrated in figures C-2 through C-4 (ref. 26). 

Tilted Dipole Model 


The tilted dipole model uses the first three terms of the spherical harmonic expansion to 
obtain an approximation of the geomagnetic field. The geomagnetic potential is 
approximated by: 

2 


To calculate the components of B at some point in space, the following relations are 
used:* 



g° P (cos 0) + (g cos X + h sin X) P. (cos 0) 

i i i * 


(C-l 9) 


B 


X 



-g° sin d + g' cos d cos X + h 1 cos d sin X 


♦Refer to appendix A for definition of geocentric coordinate system. 


(C-20) 
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(C-2 1 ) 




R„ 


\ R 

B, = - 2 


g 1 sin X - h cos X 


M 3 

R 


l cos d + g* sin 6 cosX + h 1 sin 0 sin X 


(C-22) 


The accuracy of the tilted dipole model compared to the spherical harmonic expansion 
model is illustrated in figures C-2, C-3, C-4, C-7,and C-8 (ref. 26). 
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TILTED DIPOLE 


.•* QUADRUPOLE 


EAST LONGITUDE (degrees) 


Figure C-2. Deviation of approximation models from spherical harmonic expansion model in 
equatorial plane at one R e (Earth’s surface) (ref. 26). 
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Figure C-5. -Maximum and minimum magnitude ratios of quadrupole model to the spherical harmonic 
expansion model at selected geocentric distances (ref. 26). 
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ANG 



Figure C-6. Maximum and minimum angular deviation of quadrupole model from spherical harmonic 
expansion model at selected geocentric distances (ref. 26). 
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MAGNITUDE RATIO MAGNITUDE RATIO 



SOUTH 90 60 30 0 30 60 90 NORTH 

LATITUDE (degrees) 


Figure C-7.— Maximum and minimum magnitude ratios of the tilted-dipole model to the spherical 
harmonic expansion model at selected geocentric distances (ref. 26). 
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SP-8010 (Environment) 

SP-8011 (Environment) 
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SP-8018 (Guidance and 
Control) 
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Buffeting During Launch and Exit, May 1964 

Flight-Loads Measurements During Launch and Exit, December 

1964 

Flutter, Buzz, and Divergence, July 1964 

Panel Flutter, May 1965 

Solar Electromagnetic Radiation, June 1965 

Local Steady Aerodynamic Loads During Launch and Exit, May 

1965 

Buckling of Thin-Walled Circular Cylinders, September 1965 
(Revised, August 1968) 

Prelaunch Ground Wind Loads, November 1965 
Propellent Slosh Loads, August 1968 
Models of Mars Atmosphere (1967), May 1968 
Models of Venus Atmosphere (1968), December 1968 
Natural Vibration Model Analyses, September 1968 
Entry Thermal Protection, August 1968 

Guidance and Navigation for Entry Vehicles, November 1968 

Spacecraft Magnetic Torques, March 1969 

Buckling of Thin-Walled Truncated Cones, September 1968 
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